Many-body treatment of the collisional frequency shift in fermionic atoms 
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Recent clock experiments have measured density-dependent frequency shifts in polarized fermionic alkaline- 
earth atoms using So- 3 Po Rabi spectroscopy. Here we provide a first-principles non-equilibrium theoretical 
description of the interaction frequency shifts starting from the microscopic many-body Hamiltonian. Our 
formalism describes the dependence of the frequency shift on excitation inhomogeneity, interactions, and many- 
body dynamics, provides a fundamental understanding of the effects of the measurement process, and explains 
the observed density shift data. We also propose a method to measure the second of the two 1 So- 3 Po scattering 
lengths, whose knowledge is essential for quantum information processing and quantum simulation applications. 
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Experimental efforts in cooling, trapping, and manipulat- 
ing alkaline-earth-like atoms such as Sr and Yb have led to 
unprecedented developments in optical clocks based on the 

1 S - 3 P transition G]|2lHH|5]|6l|7]|8]|9] WhUe ±C 
interrogation of a large number of atoms enhances the sensi- 
tivity of these clocks, the accompanying interatomic interac- 
tions degrade clock precision. Thus the understanding of these 
interactions is crucial for precision spectroscopy. Fermionic 
alkaline-earth atoms have also started to attract considerable 
theoretical attention in the context of quantum information 
processing ifTTl H"2l IT3l [14 1 and quantum simulation fl5l . ap- 
plications requiring the knowledge of both 1 So- 3 Po scattering 
lengths. Here we present a many-body formulation, which, 
on the one hand, allows us to understand the nature of the 
collisional frequency shift (CFS) measured in a recent exper- 
iment [3], and on the other, allows us to propose a way to 
measure the remaining 1 S'o- 3 Po scattering length, which was 
not probed in Ref. 0. Our model helps to clarify the role 
of excitation inhomogeneities, dynamics, and interactions in 
fermionic clock experiments 151 [151 [1711151 . 

Clock experiments based on Rabi interrogation start with 
a nuclear-spin-polarized sample of atoms prepared (for con- 
sistency with Ref. [3 1, which we aim to model) in an excited 
state e, which is then transferred to the ground state g by il- 
luminating the atoms during a time 1 1 with a probe beam de- 
tuned from the atomic resonance. The CFS Su> eg is inferred 
by recording the final population in g as a function of the de- 
tuning and looking for changes in the corresponding lineshape 
due to interactions. So far most treatments of CFSs in dilute 
polarized fermionic gases away from the unitarity limit were 
based on a static mean field analysis |[T6l [171 [191 l20l . The 

latter predicts a frequency shift duj eg = —jj £3 -(p g — p e )Gg 2 e , 

(2) 

with G ge the two atom correlation function at zero distance, 
which measures the probability that two particles are simulta- 
neously detected, a~ g the s-wave scattering length between the 
g and e atoms with mass M, and Pge the corresponding atom 
densities. Here we extend this formulation beyond mean- 
field and fully account for the many-body dynamics during 
Rabi interrogation. Our key statements are as follows, (i) 
Motion-induced excitation inhomogeneity can lead to s-wave 



CFS even in an initially polarized ensemble of fermions. (ii) 
CFS is sensitive to the time-averaged population difference 
between g and e atoms, consistently with the mean field ap- 
proximation, and in particular vanishes when this difference 
goes to zero (ir pulse), (iii) For a fixed pulse area, the CFS 
approaches zero as tf — > 0, meaning that, in order to experi- 
ence interactions, atoms should have enough time to feel the 
excitation inhomogeneity. (iv) Measurements of Suj eg done 
by locking the interrogation laser at fixed final ground state 
fraction are very sensitive to the pulse area and strength of in- 
teractions. Depending on these parameters, even the sign of 
5u! eg can be reversed. 

We begin our analysis with the Hamiltonian H describing 
cold fermionic alkaline-earth atoms illuminated by a linearly 
polarized laser beam with bare Rabi frequency fig and trapped 
in an external potential V(r) that is the same for g and e (i.e. at 
the "magic wavelength" JTJ). Assuming that the atoms are 
polarized in a state with nuclear spin projection mo, we omit 
the nuclear spin label and, setting H = 1, obtain |[T4l[T5l 

H = J2 Jd^l (~^V 2 + V{v)\ * Q + u- g Jd 3 rp ePg 
+lu [d 3 r(p e -p a )~Y Jd^le-^-^g + h.c.).(l) 

Here ^ a (r) is a fermionic field operator at position r for 
atoms in electronic state a — g CSq) or e ( 3 Po), while 
/5«(r) = W| t (r)W Q (r) is the corresponding density opera- 
tor. Since polarized fermions are in a symmetric nuclear state, 
their s-wave interactions are characterized by only one scatter- 
ing length a~ g , with the corresponding interaction parameter 
u~ g = AttH 2 a~ g I 'M , describing collisions between two atoms 
in the antisymmetric electronic state |— ) = (|ge) — \eg))/\/2. 
The laser with frequency ujl and wavevector k is detuned 
from the atom transition frequency ujo by 5 = lo l — uiq. 

As in the experiment of Ref. JJJ, we assume that most 
atoms are frozen along one (longitudinal) z-direction, leav- 
ing in the remaining transverse x-y plane an isotropic 2D har- 
monic oscillator with frequency oj x = uj y . We can then write 
*a(r) = 4> z {z)Y. v c av (t) U:c {x)(t) Uy {v), where (jf v and <j> v 



2 



are, respectively, the longitudinal and the transverse harmonic 
oscillator eigenmodes and & av creates a fermion in mode 
v = (v x , v y ) and electronic level a. Following Refs. ll3l [T8l . 
we assume that the probe is slightly misaligned from the in- 
direction: k = k z with \k x /k z \ -c 1. Defining 
thenfW = Sl e-^/ 2 L (r] 2 z )((f)„(x)\e ik * x \<l> v ,(x)), where 



?/i = ki J 2muj- ^ ^ are trie Lamb-Dicke parameters and L n 
are Laguerre polynomials [21], laser induced sideband transi- 
tions can be neglected if Q. v y^ v <C lo x . In this regime, H 
can be rewritten in the rotating frame as 



H 



a. 



-S n ev + E v fi av - ^2 -?r {c\„c ev + h.c) 



l/l I/ 2 1/3 f4 



where ^4 



1/1 1/2 1/4 



= # 



,c ai/ , 0^ = n £ i /,(^)^o(r?f)e-^+^)/ 2 , and 
E v are single-particle energies. 

Many interaction terms in Eq. Q can be ignored provided 
that the interaction is weaker than uo x and that the 2D oscilla- 
tor is slightly anharmonic as in the experiment of Ref. (3). 
Furthermore, we note that unless v\ = and 1/3 = V4, 
\f 4> 2 Vl ^i 3 \ > l/rij^il- As a result, the interaction is 
dominated by the terms where (u\, 1/3) is equal to (1/2, 1/4), 
{v^ B ,V2y,U2 X ,v iy ), {v2x,v± v ,v^V2 y ), or (1/4, u 2 ), which 
describe the exchange of modes along neither direction, along 
x, along y, and along both directions, respectively. Postponing 
the study of the terms exchanging modes along one direction 
only, the remaining terms conserve the number of particles per 
mode v. Assuming there is only one atom in each of N non- 
empty modes v = vn}, as is the case if all atoms 
are initially in the same internal state, H can then be reduced 
to a spin-1/2 model describing these modes: 



He 



-5S Z 



M~ J2 Uvv>(§ u ■ -1/4). (3) 



ftt 



Here U U i/' ^ e g-^uvu'u' > 5 — 2 a' Wx/-' ' ' " ' 1 ' ' 

where a are Pauli matrices in the {e,g} basis, Si =XiViZ = 
J2v 5f , an d constant terms were dropped. Hg is reminiscent 
of solid-state spin Hamiltonians, which can also feature long- 
range interactions and rich nonequilibrium dynamics l22l . 

The rotational invariance of the interaction term in Hs 
(oc U vv ' ) is key to understanding some of the basic features of 
the model. The interaction term is diagonal in the collective 
angular momentum basis \S, M, q), satisfying 5 2 |5, M, q) = 
5(5 + l)\S,M,q) and S z \S,M,q) = M\S,M,q), with 
5 = 0,... N/2 and -5 < M < 5. Here the extra label q 
is required to uniquely specify each state. The fully symmet- 
ric (Dicke) 5 = N/2 states do not interact. They are unique 
and the label q can be omitted for them. 

For a homogeneous excitation, f2„ = ti, the term 
Tl v Qv x £>x commutes with 5 2 , the interaction energy is con- 
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FIG. 1: (color online) CFS, 8uj eg / (2n), at resonant transfer, for dif- 
ferent a~ g (in Bohr radii, ao). The time-averaged ground state frac- 
tion was varied by changing t / . The solid lines were calculated by 
thermally averaging N g computed from Eq. {3}, and the dots by 
a single realization of v randomly chosen out of those satisfying 
AQ(£>) = {Att} T , = (H)t, and U(v) = (U) T . The dashed 
lines show Eq. j5} evaluated at thermally averaged parameters. Here 
N = 7,T= l^tK, (An} T /(Ct}T = 0.2 and Q = 6tt/(80 ms). The 
inset shows J <f>i 1 (x)(j>i 2 (x)dx in arbitrary units. 



served, and no CFS will be observed provided the initial state 
is an eigenstate of the interaction or a classical mixture of 
them, consistent with Ref. ifTTl . If the system is prepared 
in the Dicke manifold, states with 5 < N/2 are never pop- 
ulated and the ground state population evolves collectively 



asN^ 0) (t,S) = N 



Q. 1 

n 2 +s 2 



sin 



tVra 2 +<5 2 



where the super- 



script indicates a homogeneous excitation. 

In the presence of excitation inhomogeneity, there are 
two simple limiting cases where CFSs are absent: the non- 
interacting regime and the strongly interacting regime where 
interactions dominate over Rabi frequency inhomogeneity. 
The suppression of CFSs in the latter case is a consequence 
of the large energy gap between states with different 5, which 
brings 5-changing transitions out of resonance 11231 . In the 
intermediate interaction regime, on the contrary, the excita- 
tion inhomogeneity cannot be ignored. It will transfer atoms 
between states with different 5 generating a net CFS even at 
zero temperature. 

In most experiments, e.g. Ref. (3J, 5uj eg is measured by 
first locking the spectroscopy laser at two points, Si 2, of 
equal height in the transition lineshape (equal final ground 
state fraction under the initial condition of all atoms in 
state e) and then determining the change in the mean fre- 
quency as the interaction parameters or density are varied, 
5u! eg = (Si + 82) /2. Note that for a homogeneous excitation 
N^ 0) (t,S) = N^ 0) (t,-8) and therefore 5uj eg = 0. Defining 
Cl(V) to be the mean Rabi frequency over modes v and treat- 
ing Y^vfflv* — O ( i?) ) jS'^' as a perturbation, we write N g (tf,S) 
as Ng°\tf,S) + Ng 2 \tf,S) (the first order term vanishes), 
Taylor expand it around ±fl[ , and obtain 



Sujp 



N^(t f ,-S^)-N^(t f ,S^) 



d5 I 



(4) 
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To proceed further, we note that U vv i is a slowly varying 
function of |z/j — v[\ (i = x,y) [see Fig. [T| (inset)], except 
within a narrow range near v\ — v\. Provided ksT > Nuj z 
(which was satisfied in Ref. (3J), the occupied modes V are 
sufficiently sparse for the behavior of U ViV . to be dominated 
by its slowly varying part. Therefore, we can approximate 
U UiUj — > U(P). Under this approximation, the states with 
S = N/2 — 1, so called spin-wave states, are separated in 
energy from the Dicke states by U(p)N and are the only 
states exited to first order by the vector perturbation opera- 
tor X^(^x — &(v))S%. This allows us to obtain an analytic 
expression for 5ui eg that depends on f2„ ja . only through 
and AQ,(v), the root-mean-square Rabi frequency. This ex- 
pression is particularly simple and illuminating when evalu- 
ated at resonant population transfer (d^ — > 0): 

Aft 2 NU sin(A) sin[(NU)tf] 



Suj e 



Q 2 



A 



(NU)t f 



f(Cl,t f ,NU). 



(5) 

Here A(v) = Q(v)tf is the pulse area and f(ft,tf,NU) = 

l-n/(NU)tMtf(NU)/2)cot(A/2) Th de „ endence of 

(l-(NU) 2 /n 2 )[4,ain 2 (A/2)/A 2 -sia(A)/A] i 111C UC P C11UC "CC ui 

U, fi, and Ail on v is implied. We now make a few impor- 
tant remarks: (a) In the limit tf — > (A e — * 1), Eq. <|3j 

reproduces the mean-field expression 8ui eg — ► — ^^s^ 7 oc 

-NUGf)[tf -> 0, 5 ->■ 0] E US], (b) £w eff depends on 
the time-averaged population difference, (iV 9 — iV e )t, = 

Ijfo f (N g (T) - N e (T))dT\ 5 = 0<A n= = -^M^l and ex- 
actly vanishes at „4 = 7r when (N g ) tf = (N e ) tf - (c) For a 
fixed pulse area ^4, as tf — ► 0, the frequency shift vanishes as 
Su; eg -> -(A^/) 2 A^[7cot[y4/2]/(2y4), implying that in or- 
der to experience a CFS atoms need time to feel the excitation 
inhomogeneity. (d) At finite times tfNU > 1 and NU » Cl, 



Suj 



S m[t f (NU)] 



reproducing the expected suppression in 



eg ^ (NU)* 

the strongly interacting limit. 

So far we have assumed a fixed set of populated modes, 
V. At finite temperature, expectation values need to be calcu- 
lated by averaging over all possible combinations of modes 
{€>} weighted according to their Boltzmann factor. How- 
ever, since the quantities AQ.(y), Vt{u), and U iy) are sharply 
peaked around their thermal averages, to a good approxima- 
tion, (N 9 )t and thus (8ui eg )T can be calculated by replac- 
ing Afl(u) -► (AO) T , Cl(u) -> (n> T , and U(u) -► (C/) T . 



Here (C) T 



-E(i7)/(fe B T) 



The validity of this ap- 



Ea e- E ! i, )/<'B T ) 

proximation is demonstrated in Fig. IT] which also shows that 
Eq. <(3j is in fair agreement with Eq. 

Experimentally it is hard to measure CFSs close to resonant 
population transfer due to the small signal-to-noise ratio, and 
instead the probe laser is generally locked at a finite detuning 
0. Away from = 0, we have to consider the more gen- 
eral expression given by Eq. Q, and an intuitive interpretation 
is not straightforward. The interaction induced asymmetry in 
the lineshape not only can change the sign of the frequency 
shift, but, in general, makes CFS a very sensitive function of 
(A)t, $i , and (NU)t- In particular, the solid lines in Fig. 
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FIG. 2: (color online) CFS for different pulse areas, A. The final 
ground state fraction, N g (tf)/N, was varied by changing the de- 
tuning. Here four spatial modes at the corners of a square in the 
u x -v y plane, ui = (30,43), v 2 = (101,71), 1/3 = (30,71), and 
h>4 — (101, 43), were assumed with N = 2 atoms occupying v\ and 
i/ 2 at t = 0. The solid blue (red) lines, for UNIV. = 1 (3.25) and 
tf — 7r/f2o — 7ms, were obtained using Eq. p), which allows to 
populate 1/3 and 1/4. The empty squares (dots) were obtained using 
H s with an effective Sltf* = {l /2 and tf f = ir/tlff 1 ■ 



show the CFS as a function of the ground state fraction at tf, 
computed using Eq. Q for two atoms and four spatial modes. 
For weak repulsive interactions < (NU)t ^ (^)t (blue 
lines) and pulse areas greater than tt, Suj eg > for any detun- 
ing and approaches zero only as 5^ — ► and A — > tt (con- 
sistently with remark b above). For .4 < tt, 5oj eg changes sign 
as a function of N g (tf), and the zero crossing point moves to- 
wards smaller N g (tf) with decreasing pulse area. For stronger 
interactions (NU)t > (^)t > (red lines), while at large 
N g (tf) the sign of 5ui eg also depends on whether (A)t is 
larger or smaller than tt, at small N g (tf) the magnitude of 
the CFS becomes less sensitive to pulse area variations and 
recovers the expected negative sign for repulsive interactions 
(since at t — all atoms are in e). All these conclusions are 
consistent with the measurements reported in Ref. [3] since 
~10% variations in A over the course of a day could not be 
excluded. 

We now discuss the effect of the terms whose omission in 
the derivation of Eq. Q was not justified: the terms exchang- 
ing modes along one direction only. As shown in Fig.|2]for the 
case of two atoms occupying four modes, the spin model Hs 
(dots) reproduces the result of Eq. (|2]i (lines) very well pro- 
vided that we reduce J7q in H$ by a factor of 2: = ^o/2. 
This approximate result can be derived analytically by noting 
that the V\-v^ singlet (|— )) from the spin model is replaced 
in Eq. (|2| by the symmetric linear combination of the V\-v-i 
and the 1/3-1/4 singlets. We have checked that for N = 3 and 
4 with iV 2 modes, H s with fl^ 11 = afi (1/2 < a < 1) also 
reproduces well the results of Eq. Q. 

Conjecturing the validity of Hs with at larger as 
well, we used Hs to calculate the CFS for the parameters of 
Ref. J3l- The results are shown in Fig. [3] Our model gives 
reasonable agreement within the experimental data error bars 
for a range of scattering lengths a~ g — (100 — 500) oq with 
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FIG. 3: (color online) CFS predicted by our spin model for T = 
l(3)pK and (An) T /(fi)T = 0.15 (0.42), determined from the 
Rabi flopping curves [see inset]. The blue (pink) shaded area shows 
the uncertainty in CFS for T — l(3)^iK, assuming a variation in 
pulse area of (A)t/it = 1 ± 0.1. We used % n = 47r/80ms, 
t f = Tv/Qo ff , N — 15 [i.e. p ~ Kpcm -3 ], and a~ = 200a . 
The circles (triangles) show the T = l(3)pK experimental data of 
Ref. (3) at tf = 80 ms and p ~ 10 11 cm -3 . Inset: The shaded area 
shows Rabi flopping curves calculated at {AQ)t/{Q)t = 0.15 for 
a ea — — 600ao. The dashed black and solid blue lines correspond 
to two specific scattering lengths, a~ g — and 200ao, respectively. 
The dotted purple line is for {AQ)t / {&)t = 0.05 and a^ g — 0, 
while the blue circles are experimental data points from Ref. (3). 



Qq = 4fio • The net effect of Qq is to rescale the shift 
by a factor of four |28l . which is justified by the exchanging - 
modes corrections discussed above and by the large (up to 
a factor of 5 (24 J) uncertainty in the experimental determi- 
nation of the density. We also note that in Ref. |3| (ASI)t 
was inferred by fitting Rabi oscillations with a non-interacting 
model. However, the inset of Fig. [3j shows that interactions 
modify Rabi oscillations and that the non-interacting model 
can underestimate (Afi)r by a factor as large as three. All 
these issues combined with the experimental uncertainty in 
pulse area complicate the determination of the magnitude of 
a~ g from the current data. Nevertheless, our model at least 
suggests it to be positive, a~ > 0. 

In addition to a~ g , there is another scattering length a£ g , 
which characterizes collisions between a g and an e atom. 
a+ g collisions require a symmetric electronic state and con- 
sequently an anti-symmetric nuclear spin configuration. Thus 
the interaction Hamiltonian describing g-e collisions beyond 
the polarized regime is lfl4l[T5l 

H int = u s eg [ d 3 rp e p g +u a eg J2 J d 3 r& gm & em ,* grnl $ em . 

ram' 

(6) 

Here u s e ' g = (it+ ± u~ g )/2, s (a) stands for symmetric (anti- 
symmetric), m,m' = —I, . . . ,1 label the nuclear Zeeman lev- 
els, u+ = 4ich 2 a+/M, and p a (r) = £ m & am (r)$ am (r) . 

We now propose how to estimate a+ g experimentally. The 
method is identical to that of Ref. (3j except that probe light 
polarization should be circular instead of linear. The idea is 
to Rabi interrogate an ensemble of |e, mo) atoms with a cir- 



cularly polarized probe driving the |e, too) — \g, mo + 1) tran- 
sition. If a+ g = a~ g , then u a eg = and the dynamics of the 
system are identical to the ones described above, except that 
one needs to identify {|e, Too), \g, mo + 1)} as the spin basis. 
However, if a+ g ^ a~ g , the term proportional to u a eg will pop- 
ulate \e, mo + 1) and \g, mo). This issue can be overcome by 
applying an external magnetic field. If the applied magnetic 
field satisfies Bp^Ag ^> (U a )r, with /.in the nuclear mag- 
neton, Ag the differential g-factor between e and g l25l |2r51 , 
and (U a )T the thermally averaged antisymmetric interaction, 
then the processes populating |e, too + 1) and \g, mo) will be 
energetically suppressed, and the dynamics will be identical 
to the ones in Ref. [ 3 1 with CFS proportional to u s eg . By com- 
paring the CFS between the linearly and circularly polarized 
cases, one can in principle infer a+ g . 

In summary, fermionic clocks are sensitive to the CFS in- 
duced by excitation inhomogeneities. The CFS is sensitive 
to pulse area, detuning, and interaction strengths, but if mea- 
sured close to resonant transfer, it is in qualitative agreement 
with the expected mean field expression. To improve the ex- 
perimental resolution of the CFS, better control over pulse 
area variations is required. Rabi interrogation schemes can 
also be used to estimate af g . 

Note added in proof: While writing this paper, we learned 
about a related study of interaction frequency shifts in clock 
experiments (27). 
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